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Surface polaritons at a meta-material interface are proposed as qubits. The surface-polariton fields 
are shown to have low losses, subwavelength confinement and can demonstrate very small modal 
volume. These important properties are used to demonstatre interesting applications in quantum 
information, i.e., coherent control of weak fields and large Kerr nonlinearity at the low photon level. 

PACS numbers: 42.50.Gy, 42.25.Bs, 78.20.Ci 



I. INTRODUCTION 



In quantum information, quantum bits, or qubits, 
are the fundamental entity for encoding information 
Qubits are two state quantum systems examples of such 
systems include photon vertical and horizontal polariza- 
tion states, spin states, energy levels of an atom, etc @,0|. 
The interaction of light with atoms and the ability of co- 
herent quantum control of these interactions is a basic 
problem of quantum information science. For instance, 
quantum memory requires high coupling and reversible 
dynamics in the interaction between the photons and 
atoms in order to provide high efficient storage and re- 
trieval of light. Furthermore the carriers of quantum in- 
formation, qubits, have to be processed and manipulated 
through unitary quantum gates. The operation of some 
of these gates, e.g control phase shift operation (CZ gate), 
requires media with large Kerr nonlinearity that conven- 
tional media can not provide. Coherent control of light 
field was a subject of numerous investigations and consid- 
erable progress has been achieved in the last decades to 
develop many techniques including electromagnetically 
induced transparency (EIT) |4j, coherent population os- 
cillation, and hole burning for more extensive reviews 
the reader is referred to 

Following our recent work Q, we suggest in this re- 
view article surface polaritons as qubits for information 
processing. Surface polaritons (SP) are highly confined 
electromagnetic excitations at the interface of two me- 
dia 

Strong spatial confinement of SP leads to a huge en- 
hancement of the electromagnetic field near the me- 
dia surface. In conventional media (with permeability 
fi = 1) only the TM polarized electric SP modes can ex- 
ist. Promising properties of SP field can be realized with 
the advent of artificially fabricated negative index meta- 
materials (NIMM) with both the permittivity e < and 
permeability n < [IMl]: (1) both types (TM and TE 
mode) of polarizations exist. For quantum information 
applications, supporting both modes would be highly de- 
sirable in allowing polarization qubits, (2) low losses at 
frequency range that can be well controlled using exter- 
nal parameters, (3) highly confined field amplitude. This 



is useful for enhancing the atom-field coupling which is 
important for coherent control and for increasing optical 
depth required in quantum memory, and (4) large Kerr 
nonlinearity and order of ir phase shifts are possible to 
achieve in such media. 

Below we demonstrate these properties of SP fields on 
the interface of dielectric and negative index metamate- 
rial (NIMM). First we introduce the concept of the sur- 
face polaritons and photonics and show that these sug- 
gested qubits are capable of generating many interest- 
ing applications. In particular we demonstrate control 
of slow SP modes and the possibility of nanoscale SP 
based quantum memory devices [I4j4l9|. Coherent con- 
trol of SP qubits is achieved using EIT approach which 
provides a considerable longitudinal compression of the 
SP modes and enhancement of the interaction time with 
atoms near the dielectic/NIMM interface. Double EIT ( 
DEIT) is a promising tool to obtain large Kerr nonlinear- 
ity for free propagating single photon fields [20l - [23| j . Here 
we demonstrate large Kerr nonlinearity and cross phase 
shifts for interacting weak SP fields. A giant cross-phase 
modulation between the two SPs is achieved in a low- 
loss, sub- wavelength confinement regime. A mutual tt 
phase shift between the two SP pulses is attainable for the 
fields with a mean photon number of one, thereby open- 
ing the prospect of deterministic single-photon quantum 
logic two qubit gates for quantum computing Q. 



II. LOW LOSS SURFACE POLARITON QUBITS 

A meta-material interface Fig. [TJ supports surface po- 
laritons with both transverse magnetic (TM) and trans- 
verse electric (TE) polarization modes and these polar- 
ization degrees of freedom allow new type of qubits, i.e. 
surface polariton qubits. In the following we shall explore 
some of their basic physical properties. The surface po- 
laritons (SP) can be excited at the interface (located at 
z = 0) of two media using one of the standard tech- 
niques [9]. The first medium is assumed dielectric with 
constant electric permittivity £o£i and magnetic perme- 
ability jUoMii occupies the half space z > 0, and the sec- 
ond medium is assumed to be a NIMM medium which 
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occupies the half space z < 0, with electric permittivity 
£q£2{w) and magnetic permeability /UoM2(o;)- 



z 




FIG. 1: Surface polaritons are excited at a planar interface, 
between a dielectric in the upper half-space z > and a meta- 
material occupying the lower half-spaces z < 0. Propagation 
is in the positive x-direction.The dots above the interface rep- 
resent multi-level atoms. 




FIG. 2: Dispersion curves for surface polaritons. The solid 
line is for electric TM modes, while dashed one refers to mag- 
netic TE modes, for the set of parameters lj c — 1.37x 10 16 s _1 , 
7o = 2.73 x 10 13 s _1 , 7 m = 7c/1000, cj r = 0.145oj c , and 
F — 0.6. ei = 1.85, /ii = 1, and we fix the background 
permittivity and permeability at et = 2, [it, = 2.5. 



The SP fields propagate with frequency oj along the 
x-direction parallel to the interface with complex wave 
number (for TM polarized SP fields) 



K\\ = k 




H\£l - jJL 2 £l 



e\£2- 



(1) 



The wave numbers normal to the interface are related as; 



k\E2 + k 2 e\ = 0, 



(2) 



where k'j — 



pEjUj/c 2 is the z-component of the 
wave- vector normal to interface (j = 1,2), with the in- 
dices 1 and 2 referring to the two media. With similar 
expressions for TE polarized SP modes. 

The second medium is modeled [H, [l3[ by complex di- 
electric permittivity and magnetic permeability 



£ 2 (w) = £b - 
= (J,b 



oj(uj + vy e ) ' 
Flo 2 



(3) 
(4) 



where lo c is the electron plasma frequency, which is usu- 
ally in the ultraviolet region, and 7 (7 m ) is the electric 
(magnetic) part damping rate. F and w r are constants 
that depend on the geometry of the NIMM system. Thus 
the SP wave number is complex, its real part gives the 
SP dispersions while the imaginary part gives losses due 
to the damping terms. 

Figure ([2]) shows the dispersion relations for both TM 
(solid line) and TE (dotted line) for the set of parameters 
uj c = 1.37 x 10 16 s-\ 7 C = 2.73 x lO^s" 1 , j m = w c /1000, 
uj r = 0.145cj c , and F — 0.6. s\ — 1.85, \i\ — 1, and we fix 
the background permittivity and permeability at sj, = 2, 
/ij = 2.5. The inset shows the frequency range of interest 
between 0.1 — 0.3cli c . Near the frequency 0.14u;e, there is 



2.0 r 

1.5 



^ 1.0 



0.5 



0.0^ 
0.15 



0.16 



0.17 0.18 



0.19 0.20 



FIG. 3: (Color online) Absorption loss for surface polaritons 
as a function of frequency uj/uj c for TM ( solid line ) and 
TE ( dotted line). Losses are suppressed near frequency uio = 
0.167u; e . k, w 10 3 /m for TM modes and 10 5 /m for TE case. 



a narrow gap. We shall be interested in the frequency 
0.166w o which lies above the gap region. The imaginary 
part of equation (fj}, k(uj), yields SP loss which is shown 
in Fig. for both TM and TE polarized SP ficlds.The 
results reveal dips for the frequency loq where losses are 
suppressed i.e. k(l)o) ~ 0. For later use we mention 
here that the frequency ojo = 0.16w e will correspond or 
close to the transition wavelength in the DEIT medium 
of interest i.e. Pr:YSiO which is 606 nm. We call the 
low loss surface plasmon modes as LLSP modes. Con- 
ventional surface plasmons at a metal interface support 
TM modes only, such modes are known to be lossy. In 
the presence of a meta-material interface we have seen 
that surface polaritons exist in both TM and TE polar- 
iaztion. Furthermore these modes have low losses and 
this is only possible with meta-material interface. These 
low loss surface polaritons with two polarization degrees 
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of freedom constitute surface polariton qubits. 

The electric fields E associated with these low loss 
modes are given for mode A = TM, TE\ 



E x {r,t) = E 0iX (K l] ,z)e i 



ot) 



(•5) 



where E 0i x(k\\, z) is the SP field amplitude. In the low 
loss frequency range Kn ~ fen and the SP quantisation 
will determine the SP field amplitude which gives the SP 
coupling to atomic systems that we need to consider in 
later sections. We write the quantized SP field operator 
in the plane wave expansion as 



£M) = E / dfc ll [£oA(fc||)aA(fc||)e (ifc ii x - wt) +h.c. 



(6) 

where the creation and annihilation operators obey the 
usual commutation relation 



ax (fcy ) , a\, ) = 2it8\\> 5 (k\\ - k\ ) 



(7) 



The field amplitude E$\ is determined by matching the 
Hamiltonian in a dispersive [7], but lossles medium 



d 3 r [e i(\E\ 2 )+^m\ 2 )] 



(8) 



with the quantized Hamiltonian 



f 



h f=2J dfc||Mfc||) [a + (fc|i)a(Aj|i) + aOu)^^)] . (9) 

where / = d(fu)/duj, and the (...) indicates the time 
average. Using the equations above, we obtain the field 
amplitude ( for the TM case ) in the form 

Eo(k\\ , z) = N{k\\)9{z) + i-^sj e- fciz , 

+N{k\\)6{-z) (x - i^-z) e k2Z . (10) 

where 9(z) is the well known step function. The normal- 
ization factor N(k\\) determines the SP field amplitude 
and is given by 



L, 



N(k\\) 



D+ \ V S 



I Hu(k\\) 



_ \hl + \k£ _ 

D = £ 1 1, n Cl + £2- 



fc 2 | 2 + |fc||| 2 



M 2 



- |£il 



, |fa| 



and where the field confinement is defined as 
1 



= 



Re[fej 



J' =1,2. 



(11) 

(12) 
(13) 
(14) 

(15) 



The quantity L y is the medium size along y-direction, 
and L z is the mode length along the z-axis which is a 
function of the field confinement that is determined by 
the physical properties of the two media, viz-a-viz per- 
mittivities £1,2 (w) and permeabilities iii^[ui). We em- 
phasize that these equations are general and can be ap- 
plied to any material with arbitrary set of parameters 
appropriate for surface polaritons. Moreover, appropri- 
ate choice of materials and of frequencies of SP modes 
i.e. adjusting the pairs Ei^ui) and /ii j(w) can lead 
to large confinement ( small Q ) and thus decreases L z 
hence providing considerable enhancement of the SP field 
amplitude. This property will be used to increase the in- 
teraction coupling between the LLSP fields and atomic 
ensemble. 



III. SURFACE POLARITON CONTROL 

To control SP fields we use EIT approach, and study 
the dynamics of interaction between SP probe field and 
an ensemble of a lambda type three level atoms in 
medium 1 near the interface with a NIMM medium. The 
atom has two lower levels 1 and 2, and an upper level 
3. We assume a control field acting on the atomic tran- 
sition 2-3 while a probe field with frequency close to the 
frequency of atomic transition 1-3, and the transition 1-2 
is dipole forbidden. The probe field is taken to be TM- 
polarized SP, while the control field need not be specified 
at the moment. The control and probe fields could be of 
the same polarization or they could be different. The 
total Hamiltonian of the atom-SP field is 

H = H A + H F + H int (16) 
where free atom Hamiltonian Ha is 

H A = J2 MiPia + H^a , »wi = 0, (17) 

3 

for P r J nn = \m)jj(n\. The interaction Hamiltonian is 
given in the dipole approximation as 



H int = -J2d r E( rj ), 



(18) 



where E(rj,t) is the SP field given above at the position 
Tj , and the summation above runs over all the atoms in 
the interaction volume. 

Using the Heisenberg picture of motion we derive the 
operator equations in the limit of weak SP fields; 

d 

— o(fc||,t) = - [iw(fe||) + k(w 3 i)] o(fc||,t) 

• !»' (1 lL ' ! ''V ( 19 ) 



for g = ■ Eq / h the strength of SP dipole coupling to 
atom j and summation taken over all atoms {j} close to 
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the interface. Atomic operators P^n = \ m )jj( n \ satisfy 

+ i /" dfc|| (dsi • So/^e-^^e^ii^a^ii,*), 
|pf 2 = - (i^ + 7I1) i£ + i^fo)^, (20) 

where Sl c ( r j) = (<^32 ■ Eq/K) is Rabi frequency of the 
interaction with control SP field. The index s refers to 
the probe field while c to control field. The probe and 
control wave numbers along z-direction 



k s { c = ^jk 2 (Lo s . c ) - uj SiC 2 ei(uj StC )ni(uj SiC )/c 2 , (21) 

are calculated at the resonant frequencies W31 and W32, 
respectively. The transition frequencies are related to 
the wave number ku by the dispersion relation appro- 
priate for the type of polarization. The phenomeno- 
logical decay constants 721 and 731 are determined by 
the interactions of jth atom with other field modes and 
environment, and the cavity loss is given by «(w3i) as 
discussed in previous section. Adopting the usual EIT 
approximations @ , and performing algebraic calcula- 
tions we arrive at an equation for the slowly varying SP 
probe field A(t, x) = e ~ ik U^ x f dJfe||£fc(fc||)a(fc|| , t)e lk l x , 
and using the Fourier transformation A(y, x) = 
(27r)^ 1 J dte wt A(t,x) we find for field A the expectation 
value (A(v,x)) = A(v, x): 

[ S- - i— - ] A (v, x) = - [a(v) + k(uj 31 )] A(v, x), 
ydx v J 

/oo 
di/e -i,rt+[i ^- a( '' ) -' ,(h " ,l)]a il(i/,Q), (22) 
-00 

where v = w(fcii) — W31 is the SP probe field detuning 
from the central frequency w 3 i; the complex absorption 
coefficient a{v) is give by the expression: 

a =— -(721 - w) 

UQ(W3lj 

x / / dy&z = — . 

Jo Jo \V c (r)\ 2 - (^ + i72i)(^ + ir 3 i) 

(23) 

The basic group velocity v = v (uj 31 ) = Hj-| W3i of 

the probe SP field (without interaction with three level 
medium) and the field amplitude E are calculated at 
the central frequency W31 and A J = uj 31 — W31 is the jth 
atomic detuning from the central frequency W31. 
Spectral behavior of group velocity vo is presented in 
Fig ((3]) for TM polarized SP fields. We note that group 
velocity is close to zero near the gap edges around 0.14w e . 
In this analysis we work in a frequency range where SP 
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FIG. 4: Surface polariton group velocity (in absence of inter- 
action with three level atoms) in units c, as a function ui/uj c 
for TM. 



losses are low close to ujq , so we take the central fre- 
quency W31 ss cjo = 0.167w e which corresponds to basic 
group velocity of v g — 0.6c. 

The control field, which yields £l c , can be a freely prop- 
agating mode or an SP TE or TM field. Here we restrict 
to the latter case: |r2 c (r)| 2 = |f2| 2 e~ 2fe i z with k\ the con- 
trol field wave number in the z-direction for medium 1. 

As seen in Eq ([23]) the complex absorption coeffi- 
cient a(v) is modulated by the spatial exponential fac- 
tors so that the basic spectral properties of the SP-SP 
control are sensitive to the relative spatial behavior of 
the modes. In Eq (j23|) we have taken into account the 
inhomogeneous broadening of the resonant line by as- 
suming the weight function h(A) — (A w /n)/(A 2 + A 2 ) 
(where A w is an inhomogeneous broadening width, and 

A « A'');r 3 i = A w + 731 = , ?73l \ 7*? = is the 

spontaneous decay rate in free space, and i] is a factor 
that takes into account additional contribution to the 
line width due to the other processes e.g. inhomoge- 
neous broadening, collision etc. We assume for the opti- 
cal transitions the following typical values 73^ w 10 3 s _1 
and n « 10 6 — 10 7 in some laser crystals such as ruby 
at Helium temperature @. The real part of a(y) gives 
the absorption coefficient and the imaginary part gives 
the dispersions. In general the atomic density distribu- 
tion n(y, z) depends on both y and z directions, but here 
we shall consider only the situation where n(y, z) = no(z) 
i.e. depends on the distance from the interface. 

The atomic density no(z) can vary with distance z from 
the interface depending on the preparation of the atomic 
ensemble that is an important aspect for experiment. In 
particular we mention that optical pumping by the reso- 
nant laser fields can be used for preparing the atoms in 
the initial quantum state. As seen in Eq (|23|) the coupling 
of the atoms with probe and control fields are strong for 
atoms with spatial coordinates close to the interface. It 
is reasonable to use the atomic ensemble prepared within 
a small layer close to the interface with thickness z such 
that no{z) = no for < z < zq where no is constant. It 
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is important in this analysis also to investigate the opti- 
mal value for z where the SP control can be effectively 
realized. We can obtain the absorption coefficient in this 
case by integration over z in Eq (|23p in a closed form in 
terms of the hypergeometric function as follows; 



a(v,z ) = a (u> 31 )G (fef, k\, z , (i{v)) 



where 



a (w3i 



n L y 27T |ff | 2 
2k\ uo(w3i)r3i 



(24) 



(25) 



is the absorption coefficient of the atoms at W31 in the 
absence of the control field for zq > 1/fcf, and 



k\ k\ + k\ 1 



G(k s 1 ,k c 11 z ,l3) = ^-^{ 2 F 1 (l,^, . 

~)}. 

(26) 



where 2-F1 is the well known hyper geometric function, 
and /3 = (y + i7 2 i)(^ + ir 3 i)/|ft c | 2 . The function 
G{k\, k\, zq, (3) characterizes the spectral properties of 
the SP control and takes numerically maximum value 
of 1. 



IV. SPECTRAL CONTROL AND QUANTUM 
MEMORY OF SP 

In the presence of the control field Q Cl the main spec- 
tral behavior of a(v, zq) is sensitive to the spatial con- 
finement of SP fields expressed via the hyper geometric 
functions in Eq (f26|) . The general solution in Eq (|26|) 
has a rich spectral behavior and depends on the ratio 
between the control and probe fields wave numbers, but 
we consider few simple cases which could be very close 
for practical control of slow SP field in the limit of very 
large thickness of atomic layer i.e. zq 3> (2fcf) _1 , where 
it is possible to reduce the general hyper geometric solu- 
tions to more tractable forms. It is quite interesting that 
the spatial heterogeneity of the interaction determines 
interesting new spectral properties of the absorption co- 
efficient in different regimes of SP-control as we discuss 
in the following cases of interest that are analyzed at the 
resonant frequency vsO; 



ir 3 i 



+ iT 3 i P(v) - 1 



(1 



m - 



(27) 



For kf/kf = 0, this gives the usual spectral properties of 
the well-known EIT. 



< a >H ;c = 



■ :l V/ 3 (« / )arctanh[v / l//3( 



31 



(28) 



(hi)g = i;G 



ir 3 i 



-0(f)ln[l 



1 



\v=0 



-> p{v)\a[p(v)] v=0 -> 0. (29) 



This case , k{/k1 = 1, is most convenient for experimen- 
tal realization of SP control. 

av)|=^ = ^^)(l + ^)ln[l-^ y ]Uo 



Piy) k\lk\ 



(30) 



As it is seen from the above discussions that the ab- 
sorption of the probe SP field is highly suppressed for all 
cases in the spectral range around v w 0. We take the 
EIT medium as Pr:YSiO ( see next section for discussion 
), assuming the following parameters no « 10 20 cm~ 3 , 
L y pb 10A (where A = 606nm is the transition wavelength 
of our EIT medium ), z$ w 1/fcf ~ 1/xm, and n = 10 6 , 
we estimate l/ao(^3i) = 1/iin. The curves in Fig ([5]) 
and Fig([B|) show the absorption and dispersion profiles ( 
in units of l/ao(w3i) = l" m ) f° r the cases (ii),(iii), and 
(iv) above. We note also that the cases (iii) and (iv) of 
the SP control have more sharp spectral shape due to 
the fact that the atoms at distance zo ~> (kf) -1 make a 
large contribution to the SP control at the condition of 
weak control field amplitude. The described four cases 
demonstrate an interesting possibility to control the spec- 
tral properties of window transparency for LLSP field by 
appropriate choice of the atomic thickness z D . Further- 
more we note that the spectral properties of the window 
transparency can be varied even more extensively by us- 
ing another spatial shapes of the atomic concentration 
along z— direction. 

Spectral properties of the group velocity play an im- 
portant role in SP field control. Using Eq (|22p we find 
the group velocity of the probe SP field near w(fc||) « OJ31 



Uo(W3x) 



1 - a i;o(w3i)Im[^ 7 G'(^,fc5 ; ,2o,^/3)] 



(31) 



where the derivative -§^G is expressed through the 
same hyper geometric function. The group ve- 
locity behavior depends on the physical parameters 
ao, w 3 i, ki/kf, zo, v, O c and because of lack of space we 
only discuss the most important properties for control 
of SP field propagation. Spectral behavior of the group 
velocity v g (v) as in Eq (|31|) close to the exact two- 
photon resonance is shown in Fig ([7]) that demonstrates 
large enough reduction from ^0(^31) rs 0.6c (in the ab- 
sence of interaction with three level atom ensemble) to 
v g {v k 0) k 500m/s (due to interaction with three level 
medaium) at moderate parameters of the atomic systems 
(no = 10 20 cm -3 ) and large inhomogeneous broadening 

r 3 i - io^- 1 . 

For the pulse duration St — O.l^sec, the reduced 
group velocity leads to longitudinal pulse compression 
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FIG. 5: (a) The absorption coefficient ( real part of a(y, zo)) 
and (b) the dispersion ( imaginary part of a(v, zo) as functions 
of SP detuning v = u(k\\) — aj3i, T3i = 10 9 s _1 for fixed zo = 
1/fcf and different kf/ki ratios; red (dashed) at k{/kl = 1/2, 
blue (dotted) at kf/ki = 1 , and brown (solid) at kf/ki — 2. 
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FIG. 6: (a) The absorption coefficient ( real part of a{v, zo)) 
and (b) the dispersion ( imaginary part of a(v, zo) as functions 
of SP detuning v — ^(fc||) — W31, T3i = 10 9 s _1 for fixed k\ = 
k\ and different atomic layer thickness zo ; red (dashed) at 
zo — 0.5/fcf, blue (dotted) at zo = 1/fci , and brown (solid) 
at zo = 10/fef. 



^sp = f s 5f = 50/im that is 20 times smaller than the 
propagation length (L — 1/k(o;3i) w 1mm) of SP fields. 
Thus the SP-field pulse can be successfully stored in the 
long-lived atomic coherence P21 and retrieved by switch- 
ing the control field f2 c 14]. The demonstration of the 
possibility of such spectral manipulation and quantum 
memory of the weak LLSP field is important for quan- 
tum informatics with nanoscale plasmonics.The proposed 
manipulation can be useful for single photon and for more 
intensive light fields as well. As seen in Fig. 3, simulta- 
neous excitation of the LLSP modes with TM and TE 
polarization in the same frequency range close ui gives 
rise to quantum manipulation of the polarization qubits 
carried by the LLSP field. In particular one can store the 
polarization qubit in nanoscaled quantum memory by us- 
ing the above described scheme taking into account ad- 
ditional polarization properties of the atomic transitions. 
It is reasonable to find a possibility of most nontrivial ma- 
nipulation of the LLSP fields in nonlinear interactions of 
the TM and TE LLSP modes due to higher enhancement 
of the LLSP field near the dielectic/NIMM interface. Be- 
low we describe one of such nonlinear interaction for two 
temporally selected LLSP modes. 



V. NONLINEAR QUANTUM GATES WITH SP 
MODES 

In this section we show that weak LLSP fields can be 
used to generate large enhancement of Kerr nonlinear in- 
teraction between two LLSP pulses characterized by their 
different group velocities. For simplicity in the case of 
temporally separated LLSP mode we drop out possible 
manipulations of additional polarization degree of free- 
dom which however can be a serious subject of further 
investigations. The quest for large Kerr nonlinear coeffi- 
cient and large phase shifts has been vigorously pursued 
by many authors [2fjM3~i| over the past years and promis- 
ing success has been reported. One of the primary inter- 
est in large Kerr nonlinearity in quantum information, 
is that large nonlinearity enables the implementation of 
some quantum gates , e.g. controlled phase gate. Our 
scheme here applies to a gas system like rubidium atoms, 
as well as to solid state system like Pr:YSiO. We focus 
here on the a solid system for the following reasons; 1) 
EIT has been experimentally realized in this system by 
many authors [ijj US, HII ; and it could realize double EIT 
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FIG. 7: The SP group velocity due to the interaction with 
three level atom ensemle as functions of v = w(feii) — W31 (a) 
for fixed zo = 1/kf and different ratios fcf/fcf: red (dashed) at 
fcf/fc? = 1/2 , blue (dotted) at fcf/fef = 1 ,and brown (solid) at 
fcf/fei = 5, and (b) for the fixed ratio fcf /k\ — 1 and different 
atomic layer thicknesses: red (dashed) at zo = 0.5/ kf, blue 
(dotted) at zq = 1/fef, and brown (solid) at zo = 10/fcf. 




FIG. 8: Frequencies and fields of two slow LLSP fields a, b, of 
the control LLSP field fi c and energy diagram of rare earth 
Pr ions in YSiO crystal (Pr:YSiO). 



race the first LLSP at the medium output as it is depicted 
in Fi^TJ Energy diagram of the 5LA system and level 
scheme of Pr:YSiO are shown in Fig. [SJ As in previous 
works [28|, H^, [13] we derive nonlinear coupled equations 
for two slowly propagating LLSP fields by taking into 
account spatial confinement of interaction with resonant 
atomic systems in the spirit of 0], and we get the Kerr 
nonlinear coefficient of field b due to cross phase modu- 
lation of field a as: 



d 
Vbdt 



+ A ) E b (t, x) = ixJa{t, x)E h {t, x), (32) 



effects [24j, |2jJ |30( that are the mechanism we adopt in 
this paper to achieve large phase shifts, 2) this system is 
useful for information storage and retrieval and storage 
time of few seconds has already been experimentally re- 
ported 3) being solid state this system could avoid 
many problems that are present in gaseous systems e.g 
Doppler broadening due to the fact that atoms are locked 
into the solid with limited movement, and finally 4) the 
wavelength of this system (like the gas system) is com- 
mensurate with the current NIMM technology that has 
reached wavelength of 580 nm |34j |. 

We focus our attention only on one interesting scheme 
of the nonlinear interaction [2^, [3(| which uses 5- level 
atomic scheme (5LA) in Rb 87 that leads to uniform [35[ 
nonlinear phase shift for two slowly propagating interact- 
ing light pulses at the condition of double EIT effect. We 
adopt this 5LA scheme for our SP fields in the solid sys- 
tem Pr:YSiO. We investigate the possibility to enhance 
the nonlinear interaction between weak light fields by cre- 
ating slowly co-propagating LLSP fields characterized by 
large transverse and longitudinal confinement and by in- 
creasing the interaction time due to the EIT condition. 

Let us assume that the two interacting LLSP pulses 
are excited one by one at the interface input with slightly 
different adjusted group velocities v a ,b- Let the second 
LLSP pulse have larger group velocity v\, > v a and out- 



Xa 



k c )z ] 



\d2iEb\ 2 \di 5 E a \ 2 ) 



where I a (t, x) — \E a (t, x)\ 2 for field a, and k s a 



(33) 
is the real 



part of SP probe field (a or b) wave number along z and k c 
is the control field wave number along z. vi — w^ /(l+/9j) 
is the group velocity of i-th slowly propagating LLSP 
pulse, v^o is the group velocity of l-th LLSP pulse in 
absence of the resonant atoms, /3b = 2iTnbZo^[(k^ — 
k c )z }(\d 35 E b \ 2 )/(h 2 \n c \ 2 ), = e-y S inh(y)/y, z is 

spatial thickness of the atomic medium along z-direction, 
Q c is the Rabi frequency of the control field, rib is atomic 
density on the 3-th level, A is the spectral detuning, 
^24, d\5 are atomic dipole moments of atomic transitions 
presented in Fig. [51 E\ {I = a, b) is the electric field of l- 
th LLSP pulse which is given previously, and (...) denotes 
here averaging over the orientation of the atomic dipole 
moments. The self phase modulation (SPM) of field b, 
Xbi is ignored in the above equation since this term can be 
compensated for interferometrically as discussed in [3fj| . 
and references cited therein. The corresponding phase 
shift experienced by single quantum SP field b due to its 
nonlinear interaction with other single quantum LLSP 
pulse a as it traverses a medium of length L is given (see 
also [H|-[37j]) as <E>xpm — Xa(^s) x L/(6tv a fi) where St is 
the SP pulse temporal duration, and u> s is the Pr:YSiO 
transition frequency of interest. 
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0.20 



FIG. 9: The phase shift due to SP cross phase modulation 
in double EIT scheme (in units tt) as a function of probe 
field frequency u/uj a .The frequency of interest is around oj ~ 
0.167w e - 

For the 5LA in PrrYSiO we assume ideal EIT condi- 
tions which take place for small enough thickness of the 
atomic layer zq ps ~ l/k c « 1/zra, the transition 

wavelenght is 606nm [T^, [32|, [33[ the linewidth is about 
MHz and the detuning A = 4.5 MHz, Rabi frequency of 
control field il c — 1.5 MHz. We estimate the dipole mo- 
ments for such transition to be of the order 10 _2 eao where 
e is the electronic charge and ao is the Bohr radius. We 
have chosen the media parameters such that the transi- 
tion wavelength 606nm corresponds to SP frequency res- 
onant with PrYSio transition frequency u) s = 0.167w o 
which is close to the frequency where SP fields exhibit 
low losses and large confinement. The atomic density of 
levels is taken to be 10 20 cm~ 3 that is close to typical solid 
medium, and the medium size is assumed L w 0.3mm and 
the SP pulse temporal duration is assumed of the order 
ss 0.1/Lts. For these set of parameters which approximate 
the experimental data [19|, [32], l33j we show in Fig. [5] the 
phase shift of field b due to cross phase modulation with 
field a. The results demonstrate clearly that it is possi- 
ble to achieve all the requirements of: low losses, large 
confinement, large Kerr nonlinear coefficient and cross 
phase shift of order tt. Finally we note that with the 
nonlinear tt phase shift, the LLSP pulses can be easily 
separated at the output from each other in real experi- 
ment if the input two LLSP pulses are characterized by 
different polarizations (i.e., for TM and TE LLSP pulses). 



Whereas the analysis presented here pertains to the solid 
system Pr:YSiO, the low loss SP and their high confine- 
ment together with double EIT scheme is also applicable 
to gaseous system like Rb 87 with appropriate choice of 
parameters [3q . 



VI. CONCLUSIONS 

We have studied electric and magnetic surface polari- 
tons in half space metamaterial in contact with a di- 
electric. We have demonstrated the possibility of slow 
low loss SP fields with small longitudinal spatial size of 
the SP pulse and highly reduced group velocity. Thus 
the proposed scheme has a large potential for the quan- 
tum control of the SP-fields, and applications in com- 
pact quantum memory devices [l6l-[l3| localized near the 
surface of mate-materials. We have discussed only few 
possibilities of the SP field control that indicate the rich 
physical parameters that can be used in order to real- 
ize the most convenient conditions that suit particular 
applications of this scheme. In particular it would be in- 
teresting to analyze promising possibilities determined by 
the other types of geometry of the interface, as well as by 
using another spectral properties of the atomic systems 
and interaction with varied intensive control laser fields. 
Also we anticipate to find interesting possibilities for the 
SP field control for the atomic systems characterized by 
another spatial densities, in particular with periodic den- 
sity modulation along the SP field propagation. 

As we have found above, combining SP fields at NIMM 
interface with the DEIT mechanism yields the trifecta 
for large cross-phase modulation . These three sought- 
after properties are low loss, high field confinement, and 
large Kerr coefficients. The goal is to reach mutual phase 
shifts of tt at the single-photon level, which would have 
profound implications for quantum information technol- 
ogy 
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